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■ Abstract. Investigation of the final hadronic state properties of ultra- 
relativistics pp and Au+Au coUisions supplies information on freeze-out con- 

. ditions at RHIC and possible insights into early stages of these collisions. A 

, variety of particle spectra measured by STAR are studied within the framework 

■ of chemical and local kinetic equlibrium models. Here we present the extracted 
\ chemical and final kinetic freeze-out temperatures, strangeness saturation fac- 

QO ■ tor, final collective flow velocity, and the inferred flow velocity at chemical 

freeze-out. In light of those measurements we disscuss dynamical evoution of 
the collission system. 
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1. Introduction 



k>( \ Existence of the new state of matter - Quark Gluon Plasma (QGP) - has been 

}_j ■ predicted by Quantum Chromodynamics (QCD) long time ago. Many theoretical 

\ works were devoted to understanding of QGP properties, its possible signatures, 

and relevant phase transition conditions. Recent calculations of lattice QCD indi- 
cate a value of Tc ~ 170 MeV[T for the QGP phase transition critical temperature. 
Ultra-relativistic heavy ion collisions provide an experimental opportunity to probe 
nuclear matter under conditions which are believed to be sufficient for QGP forma- 
tion. 

The evolution of ultra-relativistic heavy-ion collisions at RHIC is generally 
thought to proceed through the following stages: 1) deconfinment of quark and glu- 
ons at high energy density in the initial stage, followed by equilibration of quarks 
and gluons; 2) phase transition and hadronization; 3) hadronic interactions and 
chemical freeze-out; 4) elastic scattering and thermal (kinetic) freeze-out. The 
bulk of the created medium appears in soft hadrons, which we study in this pa- 
per in attempt to gain understanding about the dynamics of heavy ion collisions. 
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Identified particle distributions provide valuable information about freeze-out prop- 
erties of such collisions. They also may contain important traces from the early, 
possibly partonic stage of the collisions. For example, the final collective transverse 
radial flow, due to its cumulative nature, contains any collective flow generated at 
partonic stages, and particle ratios may reflect statistical nature of the hadroniza- 
tion process. Study of spectral shapes allows to extract system temperature at 
thermal freeze-out; flavor composition reveals information about strangeness and 
baryon production and allows to measure the chemical freeze-out temperature and 
chemical potentials. 

STAR experiment has measured a variety of hadron species (tt^, K^, K^, K* , 
(/), p, p. A, A, E,E,n + Cl) in pp and Au-f Au coUisions at 200 GeV [El H |1 
The rich set of data provides an unique opportunity to investigate the freeze-out 
properties and dynamics of heavy ion collisions. 

2. Data Analysis 

Charged particles are detected in the STAR Time Projection Chamber (TPC) [0. 
TPC is surrounded by a solenoidal magnet, which provided a uniform magnetic field 
of 0.5T. Zero degree calorimeters and beam-beam counters [El provide a minimum 
bias trigger for Au-|-Au and pp collisions, respectively. The Au-|-Au events are 
divided into 9 centrality classes based on measured charged particle multiplicity 
within pseudo-rapidity jr^j < 0.5. These classes consist, from central to peripheral, 
of 0-5%, 5-10%, 10-20%, 20-30%, 30-40%, 40-50%, 50-60%, 60-70%, and 70-80% of 
the geometrical cross-section. For rear particles some of the neighboring centrality 
bins were combined to improve statistics. 

Particle identification is accomplished by measuring the ionization energy loss 
dE/dx (tt^, , p, p), decay topology together with invariant mass reconstruction 
{Kg, A, A, S, S, -I- fj), and combinatoric invariant mass reconstruction {(p, K*). 
Corrections were applied to account for tracking inefficiency, detector acceptance, 
hadronic interactions, and particle decays. Those corrections were obtained from 
embedding Monte Carlo (MC) tracks into real events at the raw data level and sub- 
sequently reconstructing these events. The propagation of single tracks is calculated 
using GEANT, a detailed description of the STAR geometry, and a realistic simu- 
lation of the TPC response [|H1 El Only pion spectra were further corrected 
for weak decay products, while the spectra of other particles represent inclusive 
measurements. 

3. Results 

3.1. Kinetic freeze- out properties 

At RHIC energies low transverse momentum spectra show clear deviation from the 
widely used rriT-exponential form. Mass-dependent hardening of particle spectra, 
characteristic of collective flow, was observed and found to increase with collision 



Freeze-out dynamics at RHIC 



3 



centrality. We use blast-wave model [ in order to describe the spectra and 
extrapolate our data outside the measured momentum range. The model provides 
hydrodynamically motivated description of the phase-space density at kinetic freeze- 
out, assuming boost-invariant longitudinal flow. It has been successfully used to 
describe pT-spectra with three parameters - kinetic freeze-out temperature Tkm, 
radial flow velocity and the exponent in the flow velocity profile (■^)". 

As shown in Figure 1 (black solid lines) this parameterization is able to simulta- 
neously describe spectra of common particle (pions, kaons and (anti)protons) with 
a single set of freeze-out parameters at each centrality. Centrality dependence of 
the extracted fit results is shown in Figure 2: the more central the collision, the 
lower extracted temperature and the higher collective radial flow velocity, reaching 
the values of Tun = 89 ± 10 MeV and (/?} = 0.59 ± 0.05 in the 5% most central 
coUisions [EJ. 




Fig. 1. Identified particle spectra in central Au-t-Au events Measured by the STAR 
experiment. Black solid lines show blast-wave fit results for charged pions, kaons 
and protons. Dashed lines represent blast-wave prediction for multi-strange baryon 
spectra shapes based on freeze-out parameters of tt, K, p. Red lines are blast-wave 
fits to the S and spectra. 

Blast-wave model can also reproduce spectral shapes for other particles mea- 
sured by STAR, but it appears that different set of freeze-out parameters is needed 
to describe rare particle distributions. Dashed lines shown in Figure 1 illustrate 
that blast-wave model fit to the common particle spectra fails to reproduce spectra 
shape of the rare particles. Red lines show the blast-wave model fir to the S and O 
spectra, yielding a set of parameters different from that of the common particles. 
The fit parameters of such single-spectra fits are shown in Figure 2 in colored sym- 



4 



O. Barannikova 



bols. Possible correlation of the lower values of radial flow velocity for rare particles 
with their low hadronic interaction cross-sections suggests sequential freeze-out of 
particle species, with multi-strange baryons decoupling from the system earlier than 
the common particles, that are more prone to elastic rescattering. 

3.2. Chemical freeze- out properties 

To study flavor composition at chemical freeze-out we extract integrated particle 
yields extrapolating our measurement over the whole phases-pace. In the framework 
of a chemical-equilibrium model [ 1121 integrated yield ratios can be described 
by a set of parameters: the chemical freeze-out temperature (Teh), the baryon and 
strangeness chemical potentials {^b, IJ-s), and the strangeness suppression factor 
(7s). We fit our ratios with such a model to extract these parameters. The value 
obtained for the chemical potential is independent of centrality within errors, and 
fis is consistent with 0. The obtained 7^ increases quickly from peripheral to mid- 
central collisions approaching unity for most central events. In the framework of 
the model used this is an indication of s-quark equilibration achieved in central 
heavy- ion collisions at RHIC. 

In order to gain insights into collision evolution, we contrast Tch in Fig. 2 
with the fitted kinetic freeze-out temperatures. Contrary to Tkim no centrality 
dependence is observed for Tch- Moreover, the extracted Tch value is close to the 
QCD prediction of critical temperature. As hadronization is a universal process, 
it is suggestive that chemical freeze-out in heavy ion collisions at RHIC may occur 
very close to hadronization [ I14j. In addition, multi-strange baryons appear to 
decouple earlier from the system, possibly at/right after chemical freeze-out. This 
allows to estimate radial fiow velocity at chemical freeze-out/hadronization to be 
on the order of 0.45 ± 0.10 c. Such a significant amount of radial flow has to come 
from early partonic stage of the collision, if indeed present at hadronization. 

3.3. Inferred initial conditions 

To characterize initial condition we construct the following variable: (''■^/^y^^ ^ where 
S is an estimate of the transverse overlap area based on the number of participants, 
calculated by the MC Glauber model [EI- The '^'^^/^y^^ jg related to the initial 
conditions of the collision such as energy density [[T5|. It is also a relevant quan- 
tity in parton saturation picture, which suggests that in high energy collisions the 
initial gluon density is saturated up to a momentum scale that is proportional to 

\J ^'^'^/^y 'h. [ 16 . The centrality dependence of \J ('^^A'^ lZ jg ghown in Figure 2, 
where one sees a strong increase as function of collision centrality. In other words, 
higher energy densities are achieved in more central collisions. However, despite 
of the significant change in the initial conditions, the collision system yet always 
evolve towards the same chemical freeze-out temperature. 

Using measured pion, kaon and (anti)proton spectra we estimate a mean trans- 
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Fig. 2 
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(a) J i^mi^mliL (stars), fch (circles) and Tkm (triangles); (b) (/?) and 

(c) 7s as a function of the charged hadron multiplicity. Black symbols represent 
results extracted from tt, K, p measurements and shown with systematic errors. 
Preliminary results for rare particles are shown in colored symbols (as marked on 
the plot), error bars are statistical only. 



verse energy per particle in the 5% most central events to be (Et) ~ 600 MeV. At 
RHIC energies the pion multiplicity in the final state correspond closely to the initial 
gluon multiplicity [El- This would imply an initial temperature of T « 300 MeV 
assuming the ideal relativistic gas dynamics. For the case of 2 quark flavors, such 
temperature would give an initial energy density e = ^^|^T^ 10 GeV This is 
in agreement with Bjorken energy density estimates by STAR [^j and PHENIX [ 
[T^ . and is significantly above the QCD predicted critical energy density, support- 
ing that a new state of matter created at RHIC in the early stage of the Au+Au 
collisions. 



4. Conclusions 

The Study of a variety of hadron spectra from STAR within the framework of chem- 
ical and local kinetic equilibrium model yields the following conclusions: 1) initial 
energy density increases with centrality and is sufficiently above the QCD predicted 
critical energy density. 2) Chemical freeze-out seems to occur at a universal temper- 
ature, close to the predicted critical temperature for the phase transition. 3) Rare 
particles appear to freeze-out kinetically with significant amount of flow at/right 
after chemical freeze-out. 4) The bulk of the medium evolves further towards final 
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kinetic freeze-out with lower temperature and larger flow velocity in more central 
collisions. 
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